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Abstract: 

There are two types of plume flow models: (i) 1 D models using ad hoc spreading 
functions, f(r). (ii) MHD models. 1 D models can be multifluid, time dependent, 
and incorporate very general descriptions of the energetics. They confirm 
empirical results that plume flow is slow relative to requirements for high speed 
wind. But, no published 1 D model incorporates the rapid local spreading at the 
base (fl(r)) which has an important effect on mass flux. The one published MHD 
model is isothermal, but confirms that if b=8pp/|B|2«1 then the field is nearly 
potential below -70,000 km. Building on the MHD result, we apply a two scale 
approximation to calculate fl(r). We also compute the global spreading (fg(r)) out 
to 5.0 RSUN imposed by coronal hole geometry. Global MHD models provide a 
potent method of calculating fg(r). Unambiguous plume signatures have not yet 
been found in the solar wind. This is probably due to strong mixing of plume and 
interplume flows near the Sun. We describe a physical source for strong mixing 
due to the observed flows being unstable to shear instabilities that lead to rapid 
disruption. 



DRAFT, 18 February 1998 


MODELS OF PLUMES: THEIR FLOW, THEIR GEOMETRIC SPREADING, AND THEIR 

MIXING WITH INTERPLUME FLOW 


S. T. Sues s 

NASA Marshall Space Flight Center/ES82, Huntvsille, AL 35812, USA 


ABSTRACT 

There are two types of plume flow models: (i) ID 
models using ad hoc spreading functions, f(r). (ii) 
MHD models. ID models can be multifluid, time 
dependent, and incorporate very general descriptions cf 
the energetics. They confirm empirical results that 
plume flow is slow relative to requirements for high 
speed wind. But, no published ID model incorporates 
the rapid local spreading at the base (f(r)) which has an 
important effect on mass flux. The one published MHD 
model is isothermal, but confirms that if 
P=8rcp/|B| 2 «1 then the field is nearly potential below 
-70,000 km. Building on the MHD result, we apply a 
two scale approximation to calculate f(r). We also 
compute the global spreading (f g (r )) out to 5.0 Rsun 
imposed by coronal hole geometry. Global MHD 
models provide a potent method of calculating^^). 

Unambiguous plume signatures have not yet been 
found in the solar wind. This is probably due to strong 
mixing of plume and interplume flows near the Sun. 
We describe a physical source for strong mixing due to 
the observed flows being unstable to shear instabilities 
that lead to rapid disruption. 

Key words: Plumes, Theory, MHD, Kelvin-Helmholtz 


1. INTRODUCTION 

Plumes are bright quasi-radial rays between one and 
several Rsun in coronal holes. They are observed during 
eclipses and, e.g., from SOHO and SPARTAN-201. 
Using SOHO, DeForest et al. (1998) have finally 
unambiguously shown all plumes lie over photospheric 
magnetic flux concentrations, although not all flux 
concentrations have plumes. Modeling, SOHOAJVCS, 
SPARTAN-201, and Interplanetary Scintillations 
(IPS) have shown that plume flow speeds are typically 
-200-300 km/s at 2.0-5.0 R S un (Wang, 1994, Habbal 
et al., 1995; Corti et al., 1997; Poletto et al., 1997) 
while the interplume flow speed may already be 750 
km/s at 5.5 Rsun (Grail et al., 1996). Therefore, plumes 
flow much more slowly than interplume plasma inside 
10 Rsun. 

Being bright in white light, plumes are denser than the 
interplume plasma. This, together with the information 


given above, suggests that plumes should be 
observable in the interplanetary medium. However, 
concerted searches of the Ulysses data (McComas et al., 
1995; Neugebauer et al., 1995) have yielded 
inconclusive results. In fact, the high speed solar wind 
coming from coronal holes is remarkably smooth 
(Phillips et al., 1995). Therefore, plume and 
interplume plasma are mixed somewhere close to the 
Sun. 



Figure 1. Plumes lie over some , not all, photospheric 
magnetic flux concentrations. Spreading is therefore 
rapid near the base and governed by the coronal hole 
geometry above - 50,000 km (see exaggerated plume 
geometry at right). 

Modeling plumes therefore reduces to at least three 
individual problems: the flow, the geometric spreading, 
and the plume/interplume mixing. The flow modeling 
is fundamentally straightforward and has a good 
foundation. The directions in which it should go are 
clear. It is easy to model the field, or streamline, 
geometry of plumes because P«1 in coronal holes, out 
to at least 10 Rsun and probably to 0.5 AU. Given the 
streamline geometry, one-dimensional (ID) solar wind 
models can be used to determine the flow properties in 
plumes using spreading functions. However, the 
plume/interplume mixing problem must be resolved. 
Plumes can be observed with LASCO to at least 10 
Rsun, but are well-mixed by 0.3 AU. 

In this brief summary, flow models are reviewed first, 
then the technique for calculating the streamline or field 
geometry is discussed - since parts of it are quite new, 
and finally an idea is outlined for plume/interplume 
mixing in coronal holes to produce what is observed 
farther from the Sun. 
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2. PLUME FLOW MODELS 

For the purpose of calculating plume flow in the 
corona, inside ~10 Rsun, an acceptable assumption is 
that the flow is radial and that the plume geometry is 
expressed in terms of the “spreading function” 
introduced by Kopp and Holzer (1977). Section 3 will 
address calculating the spreading function. With these 
assumptions, the flow calculation reduces to a ID 
problem. Specifically, this leaves only the momentum, 
continuity, and energy equations along the streamline, 
with multiple momentum and/or energy equations fix 
multicomponent flow, plus the equation(s) of state. 

In spherical coordinates, where A a is the area at the base 
of a streamline or flux tube, at radius r a , f(r) is defined 
in terms of A(r) by 

A(r) = (r/r 0 ) 2 f(r)A o (1) 

Many methods are used to solve the ID, possibly time 
dependent solar wind equations. A particularly good 
method is as an “initial-boundary value problem” in 
which an essentially arbitrary initial state is allowed to 
relax in time until the solution is steady. This 
correctly deals with the critical point(s), a recovered 
solution is guaranteed to be stable, and it can be very 
efficient, numerically, if implicit numerical time 
differencing is used. A time dependent solution has the 
added benefit that it can also be used for simulations of 
transients. Various realizations of this approach are in 
the literature. Two fully implicit, noniterative models 
arc those by Suess (1982a) and Hu et al. (1997). The 
later paper poses a fairly complex problem cf 
anisotropic temperatures, heat and momentum sources, 
Alfvdn waves, an electron-proton fluid, and inhibition 
of radial thermal conduction by the Archimedian spiral 
of the interplanetary magnetic field. The calculation 
was only applied to fast wind analysis, but could 
equally well have been applied to plumes. The 
differences between interplume and plume flow lie in 
the densities chosen, the heating terms, and in the 
geometry near the base of the flow. 

A particular model is that by Habbal et al. (1995), 
discussing “less dense” and “denser” structures in 
coronal holes and which seems to apply to plumes. It 
is strongly constrained by SPARTAN 201-01 
observations made of the northern polar coronal hole in 
April 1993. The model is two fluid (electrons and 
protons) and has Alfv6n waves and radiative losses 
(Rosner et al., 1978). Radial flow is invoked in the 
plumes, while the spreading function is greater than 
unity in the interplume flow. The results are that flow 
speeds in plumes are less than in interplume flow. The 
smaller spreading function in the plumes helped 
produce the required enhanced density and was partly 
responsible for the smaller flow speed. Plots for their 
flow speeds in the interplume and plume are shown in 
Figure 2. 




Figure 2. Flow speeds from Habbal et al (1995): □ r s 
and X’s are empirical values based on observed 
density , assumed spreading function , and a mass flux 
at 1 AU of 1.6 or 3 x l(f cm 3 . The solid lines are the 
model computations. Left: ambient coronal hole . 
Right: dense structures (“plumes ”). 

Another particular model is that by Wang (1994), who 
solved the steady state equations, also produced low 
speed, dense plumes, but differed in the details of the 
physics from that of Habbal et al. Wang’s one fluid 
model used equal spreading functions in the plume and 
interplume. He constrained his model to give realistic 
velocities, temperatures, and mass fluxes at 1 AU, had 
heating and radiative losses in the corona and transition 
region, and used extra heating in plumes to produce 
their enhanced densities. 

Wang’s and Habbal et al.’s models are examples cf 
how ID solar wind solutions can be applied equally 
well to plumes and the standard solar wind. Both 
models produce speeds similar to empirical results 
(Corti et al., 1997; Poletto et al, 1997) at ~2 Rsun and 
have densities higher than interplume plasma. 
However, from this point they diverge. Habbal et al. 
use plume/interplume spreading differences to help 
produce enhanced the plume density while Wang uses 
enhanced basal heating to do the same. Neither model 
invokes the rapid geometric spreading now known to 
exist at the base of plumes (DeForest et al., 1998) and 
which has been shown to be dynamically important 
(Del Zanna et al., 1997). Finally, the interplume flow 
speed is less than the reported 750 km/s at 5.5 Rsun 
from IPS observations (Grail et al,, 1996). 


3. PLUME GEOMETRY MODELS 

The geometric spreading of plumes can be computed 
with acceptable accuracy independently of the flow 
because P«1 in plumes and throughout the 
surrounding coronal holes from the base of the corona 
to at least 10 Rsun. This has been explicitly proven by 
Del Zanna et al. (1997), and was invoked in earlier 
potential field models of Newkirk and Harvey (1968) 
and Suess (1982b). Del Zanna et al. used an isothermal 
model to show that the flow has negligible effect on 
geometric spreading. Their results are summarized in 





Figure 5. Array of flux concentrations at the photosphere , with rapidly spreading field lines above the photosphere to 
a height — the typical flux concentration separation distance . Some concentrations have plumes, but not all Above 
SO, 000 km, the field is locally smooth in the transverse direction. 
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By a height comparable to the typical separation 
between flux concentrations, the field is smooth and is 
dominated by the large scale geometry of the coronal 
hole. This may seem counterintuitive, so a brief 
explanation is given. Coronal holes are open, releasing 
solar wind, because of the radial pressure gradient. The 
field is too weak to contain the plasma, but the plasma 
P remains small until at least 1 0 Rsun because of the 
low density in coronal holes. This means that 
transverse (to the radial direction) gradients in the field 
will be rapidly smoothed but that transverse gradients 
in the plasma density have little effect on the geometry 
(Suess & Smith, 1996; Suess et al, 1998a). This is the 
reason Ulysses measured essentially constant radial 
magnetic field strength across the polar coronal holes 
and it is why plumes and adjacent interplume 
streamlines have the same geometric spreading above 
the rapid basal plume spreading region. This permits a 
very simple separation to be made to compute the 
spreading function - a two scale approximation. Let 

f( r ) ~~ flocal( r ) fgiobal(r) & flfrffgfr) (2) 

Then 

B(r) = BofRsus / rf l/(fi(r)f t (r) (3) 

With the two scale approximation, f(r) is computed 
from a potential field model incorporating the 
distribution of magnetic flux concentrations at the 
photosphere, and f g (r) is computed from a global 
model of the corona such as that illustrated in Figure 6 
(Wang et al., 1998). The calculation is described in 
detail by Suess et al. (1998a). 



Figure 6, The magnetic field lines in the MHD global 
coronal model of Wang et al (1998). 


Figure 5 shows a plot of the array of flux concentrations 
at the base of the corona - some with plumes. The 
results shown below are for this rectangular 
distribution. A hexagonal distribution might be more 
appropriate for the Sun but the quantitative difference, 
in terms of the spreading function, is negligible (Suess 
et al., 1998a). A calculation with a 5% background 
flux, full width at half-maximum of the concentration of 
about 0.12 of the distance between concentrations, 


results in the f(r) shown in Figure 7. fi(r) increases 
rapidly from unity to -13 between the photosphere and 
<50,000 km. This calculation was scaled such that the 
concentrations were separated by -30,000 km. Above 
50,000 km fi(r) is constant. This is a requirement cf 
the two-scale approximation, and is well satisfied. 


<■) (b) 



Figure 7. Combined local and global spreading 
factors for a flux concentration of half-width-0.13 
times the concentration separation distance and a 5% 
background field. The global spreading is computed 
along the symmetry axis in the Wang et al. (1998) 
model (from Suess et al, 1998a). 

For f g (r), the geometry is determined by the coronal 
hole geometry. This can either be estimated 
empirically, by examining coronal hole boundaries 
(Munro & Jackson, 1977) or theoretically, from global 
MHD models. MHD models show that the spreading 
depends only weakly on position across a hole, except 
very near the edges, so that empirical estimates have 
considerable validity. However, here we will show 
results from the MHD calculation of Wang et al. 
(1998). The Wang et al model incorporates volumetric 
source terms and latitude dependent boundary 
conditions to obtain physically realistic flow speeds 
and densities in the coronal holes and reasonable 
temperatures and densities in the streamers. It gives a 
plasma P throughout the coronal hole of the same 
magnitude as estimated by Suess and Smith (1996) 
and Del Zanna et al. (1997). The field lines are shown 
in Figure 6 for this model. We have computed f g (r) 
from the model for field lines near the edge and at the 
center of the hole. The spreading functions derived from 
this model are qualitatively similar to those derived 
from Steinolfson et al.’s (1982) model and from Suess 
et al.*s (1996) model, as well as those derived 
empirically. f g (r) along the axial field line in the model 
shown in Figure 6 is that plotted in Figure 7. It is seen 
to vary only slowly between 1.0 and 1.2 Rsun, and to 
vary much more rapidly than f(r) above 1,2 Rsun. This 
is the second requirement for the applicability of the 
two scale approximation, and it is also well satisfied. 
The product of f(r) and f g (r) gives the total spreading 
function, which is also plotted in Figure 7. Because cf 
the rapid basal spreading, the total spreading ultimately 
reaches almost 50 at 5 Rsun. The background field is 
very important in this model, with the local spreading 
reaching 50 for zero background. 
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4. MIXING OF PLUME & INTERPLUME FLOWS 

Plume and interplume flow speeds are inferred to differ 
by at least 200 km/s at 5 Rsun from empirical evidence 
(Grail et al., 1996; Corti et al., 1997; Poletto et al., 
1998) and from models (Habbal et al., 1995; Wang, 
1994). The models also suggest that the temperatures 
and densities could differ at 1 AU. These signatures 
would be easily found in Ulysses measurements cf 
polar coronal hole flow. The rotational (corotating 
interaction region) interactions which otherwise 
obliterate plume signatures in equatorial flow arc 
almost entirely absent at high latitude (Neugebauer et 
al., 1995). Yet, no obvious plume signatures have been 
detected. This simple empirical feet leads to the 
conclusion that there is mixing of plume and 
interplume plasmas somewhere between 10 Rsun and 
0.3 AU (the perihelion of Helios 1/2). 

There are probably numerous processes which could 
lead to plume/interplume mixing. It is also possible 
that momentum and energy deposition in the outer 
corona could erase plume/interplume differences. 
Nevertheless, whatever other processes exist, it is easy 
to show that plumes will be subject to MHD Kelvin- 
Helmholtz (KH) shear instabilities beginning at ~10 
Rsun, and that these instabilities will otherwise lead to 
disruption of the plumes and mixing with interplume 
plasma. It is furthermore of interest that, because this 
occurs in a low to moderate P plasma, the KH 
instability will produce Alfvdnic fluctuations - 
potentially a source for some of the AlfV6nic 
fluctuations observed in the solar wind. This is a 
strong hypothesis, based on linear stability and growth 
rate analysis and on extensive numerical simulations. 
However, no published simulations or detailed 
evaluations explicitly address parameters appropriate for 
coronal plumes so considerable analysis remains to be 
done with respect to this hypothesis. The physical 
process and existing numerical results are outlined 
here. 

Briefly, the ordinary KH instability occurs when the 
shear speed between to fluids becomes too large. It is 
reviewed by Chandrasekhar (1961), and can occur even 
in an incompressible, inviscid fluid in which there is 
horizontal streaming. In the viscous fluid case, the 
interface is unstable for ail wave numbers greater than a 
specific value determined by the velocity difference, no 
matter how small the velocity difference. For 
compressible fluids, this narrows to the shear having to 
be larger than the sound speed for the KH instability to 
occur. 

Just as one might intuitively expect, the MHD KH 
instability occurs for shears greater than the Alfven 
speed (Hardee et al., 1992; Hardee, 1995; Hardee & 
Clarke, 1995). This has been examined in detail for 
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linear stability in slab jets and in numerical 
simulations for slab jets and jets in cylindrical, 
expanding atmospheres. It has been found that the 
instability grows most rapidly for AlfV6nic (transverse, 
or sinusoidal) fluctuations and that nonlinearities do 
not stabilize the instabilities. "The MHD KH results 
suggest that a jet which is initially subAlfvdnic and 
stable to disruption will be doomed to disruption at 
the Alfven point if it becomes superAlfvdnic as a result 
of jet expansion” (Hardee et al, 1992). 

Turning to coronal holes and plumes, the first thing to 
consider is whether the stability criterion of the MHD 
KH instability has any relevance. Figure 8 (Krogulec et 
al., 1994) show one example of estimating the Alfven 
speed in coronal holes. In the right panel, it is seen 
that the AlfVdn speed peaks above 1000 km/s and is 
generally above several hundred km/s out to 10 Rsun. 
Then, for typical field strengths, the Alfven speed drops 
below -400 km/s, which is typical of the expected 
shear between plume and interplume plasma. By way 
of confirmation, the Alfv6n speed shown here is similar 
to that estimated independently by Suess (1988) in an 
independent model. 

Using Figure 8 and the above discussion of the values 
of P in coronal holes, it is possible to construct the 
cartoon in Figure 9 illustrating the parameter regimes 
which plumes typically pass through. For moderate 
shear (>,~ 300 km/s) plumes first become 

superAlfV6nic and then move into a high P plasma 
regime. For small shears, it is possible that plumes 
would move into the high P regime before becoming 
shear unstable and, hence, the resulting KH instability 
will be more like die OHD KH instability. 

The remaining important parameters governing the 
behavior of the KH instability are the internal and 
external: magnetosonic Mach numbers, Alfv6n Mach 
numbers, and densities. 

Several MHD simulations have been performed by 
Hardee and his colleagues. However, these have all 
been directed towards the analysis of astrophysical jets. 
Although they confirm the linear stability analysis, 
they do not incorporate parameters appropriate for 
plumes. In particular, none of the reported simulations 
have a magnetic field outside the jet Nevertheless, it is 
instructive to examine the results shown in Figure 10 
from one of the simulations. The simulation is for a 
slab jet in a cylindrical geometry, with an expanding 
atmosphere. The geometry is independent of the 
coordinate normal to the plane of the figure. Due to the 
expanding atmosphere, the shear, which is initially 
subAlfv£nic, eventually becomes super AlfV6nic. At this 
point, the MHD KH instability occurs and sinusoidal 
oscillation set in. The plot shows the resulting 
magnetic field lines. 
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Figure 8 . Temperature distirbution (a), electron density (b), and flow and Aljven speeds (c) for a model coronal hole 
atmosphere (Krogulec et al, 1994). The Aljven speed has been computed for several difference photospheric magnetic 
field strengths. 



Figure 9. Different plume regimes in shear velocity flvi-v e |) and f3 in the corona, v* is the Aljven speed. 7 ’ and *e* 
refer to internal and external to the plume. The plume is initially subAlfvinic in a low-fi plasma. For moderate shear 
(>,-300 km/s) it first becomes superAlfvinic and then moves into a high-f3 plasma. 



Figure 10. Magnetic field lines for a jet moving into 
an expanding external atmosphere in cylindrical 
geometry (Hardee and Clarke, 1995). The jet is 
magnetized and remains stable until it becomes 


super Aljvenic. It then destabilizes abruptly. There is 
no field in the external medium. 

Presently, simulations are under way for conditions 
more appropriate for coronal plumes (Parhi & Suess, 
1998; Suess et al., 1998b). These simulations use the 
same numerical code used by Hardee, the so-called 
ZEUS code. The present version is fully 3D, time 
dependent, and contains the metrics for calculations in 
Cartesian, cylindrical, or spherical coordinate systems. 
Preliminary simulations are being made only in 
Cartesian slab geometries in order to carefully map out 
the parameter regime. 


5. SUMMARY 

Plume flow can be well-simulated by applying standard 
ID, time dependent solar wind models. Several 
appropriate models already exist, including those of Hu 
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(1997) and Habbal et al. (1995). Plume geometry is 
incorporated through the use of a spreading function. 
The differences between plume and interplume solar 
wind models are in the heating and momentum source 
terms and in the spreading function near the base. 

The spreading function is calculated by combining a 
potential field calculation at the base, to account for the 
magnetic flux concentrations located there, with a 
global MHD coronal model above heights of -50,000 
km. Plume and adjacent interplume spreading functions 
are essentially identical between 50,000 km and where 
the plasma P becomes 0[ 1] or where the plume starts 
mixing with the interplume. The total spreading 
between the base of plumes and the interplanetary 
medium is between 30 for large coronal holes and -200 
far small coronal holes. This also depends on the 
background magnetic field in between photospheric 
magnetic flux concentrations. A 5% background (in 
strength) is used here. 

The inferred mixing of plume and interplume plasma is 
not a resolved process. However, such ideas as the 
“pulsed solar wind model” of Feldman et al. (1998) do 
not apply since they only address mixing along flux 
tubes, not across flux tubes. Plumes will be subject to 
the MHD KH instability if the shear speed is greater 
than the AIfv6n speed. Because the ambient Alfv6n 
speed decreases with increasing heliocentric distance, 
becoming generally <50 km/s at 1 AU, almost any 
shear will eventually become unstable. It appears that 
shear associated with plumes will become unstable at 
15 ±10 Rsun, generating Alfvenic fluctuations because 
the ambient P«1 and becoming a potential source for 
such fluctuations seen further from the Sun in the solar 
wind. 


6. UNRESOLVED ISSUES 

Although plume modeling is under way, it has not 

gone very far. Much more is needed. An incomplete 

list of useful additions is offered here: 

1. Application of the best ID models to plumes (i.e. 
dense, slow wind). 

2. Incorporation of rapid basal spreading, shown by 
Del Zanna et al. (1997) to have an important effect 
on mass flux. 

3 . Parameteric survey of consequences of changing the 
photospheric field morphology and strength. 

4. Time dependences: This category requires some 
comment. Many time fluctuations have been 
reliably reported at the base of plumes and, more 
anecdotally, moving up plumes. Bright points and 
flaring bright points sometimes lie at the base cf 
plumes. It seems obvious that transients in plumes 
are an open area for modeling but nothing has been 
done. Since the most user-friendly plumes models 
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are inherently time dependent, much could be 
learned by using them to study transients. 

5, Finally, plume/interplume mixing is an open area 
for study. 
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